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J£} ■ ABSTRACT 



Whatever the ultimate energy source of gamma-ray bursts turns out to be, the 
resulting sequence of physical events is likely to lead to a fairly generic, almost 
unavoidable scenario: a relativistic fireball that dissipates its energy after it has 
Ph \ become optically thin. This is expected both for cosmological and halo distances. 

Here we explore the observational motivation of this scenario, and the consequences 
Cn| \ of the resulting models for the photon production in different wavebands, the ener- 

getics and the time structure of classical gamma-ray bursters. 

> 

O ■ 1. OBSERVATIONAL REQUIREMENTS 

o\ '. 

A main requirement of almost any model for gamma-ray burst sources (GRB), 
CN | which has been long appreciated, is that the "working surface" must be expanding 

due to radiation pressure. This is expected either for an extended galactic halo 
(D ^ 10 24 cm) or for a cosmological (D ~ 10 28 cm) spatial distribution (distances 
much smaller than 100 kpc are not favored by current observational limits on the 
isotropy of the location distribution, e.g. Fishman, 1994). For a characteristic 
observed fluence F = 10 _6 -F_6 ergs cm -2 the energy in a solid angle 9 2 is 



E = 10 43 F_ 6 D% 4 9 2 = 10 51 F_ 6 D% 8 9 2 ergs, (1) 



so the luminosity is highly super-Eddington, E/tf, S> 10 38 (M/M Q ) ergs for the 
characteristic burst durations if, <> s and masses M ~ Mq expected for most 
^ ' progenitors. Coupled with the fact that gamma-ray photons above e ± pair thresh- 

old constitute most of the observed flux, this led to the concept of an expanding 
pair-photon fireball (Cavallo and Rees, 1978, Paczyhski , 1986, Goodman, 1986). 
However, from the above information alone it is not possible to determine whether 
the expansion will be slow or fast. If the baryon load of the radiation-pressure 
ejected shells were large (E/Mc 2 & 1, the expansion would be subrelativistic, and 
the flow would remain optically thick, leading to a degradation of the gamma-rays 
(Paczyhski , 1990). 

A second element which must enter any successful GRB scenario is implied by 
the detection in many GRB of photons with energies e 7 *S 1 GeV. Pair formation 
sets in for 

e 7 > 2(m e c 2 ) 2 [E t {\ - cosa)] -1 ~ 4(m e c 2 ) 2 'e^ 1 a~ 2 (2) 

in the laboratory frame, where et is the lab frame target photon energy and a is 
the relative angle of the two photons. However in the same frame, causality implies 
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a £s r 1 where V is the bulk Lorentz factor. Therefore V 1 ^ a 2m e c 2 / ^/IJe^ , 
or since e t ~ 1 MeV is where much of the lab-frame GRB photons are, 

6 7 , MeV ^ lOVMevil Me V . (3) 

This (e.g. Fenimore, Ho and Epstein, 1993, Harding and Baring, 1994) is sometimes 
referred to as the need for a "beaming" characterized by a factor V. However, it 
must be stressed that this refers only to relative photon angles - the GRB could 
perfectly well be emitting isotropically. (The word "beaming" can be misleading, 
and is better reserved for actual channeling of the total emission into a restricted 
solid angle). In any case, one infers a highly relativistic expansion, r > 1, and 
this in turn implies that, somehow, the GRB deposits much of its radiation energy 
in a low density region where it significantly exceeds the baryon rest mass energy 
density. This "low baryon loading" is required by the observation of high energy 
photons e 7 ^ 0.1 — 1 GeV. 

A third requirement of a satisfactory GRB model is that it must account for the 
generally non-thermal appearance of the spectra, strongly suggesting an optically 
thin spectrum from power-law relativistic electrons. In principle a thermal electron 
distribution in a scattering-thick medium could produce a power law up to an energy 
~ 3kT e , but in order to have the photon power law extend to GeV energies the 
"temperature" would have to be of the same order, and at these energies nonthermal 
distributions are more likely. The requirement is therefore very likely to be that the 
emitting plasma is optically thin and nonthermal. The comoving electron density 
is n' e = L/4irr 2 c 3 rir, where primes denote comoving- frame (CF) quantities, and the 
scattering depth r e ~ n' e aTr /V must satisfy 

r e ~ La T /(4irrm p c 3 r]T 2 ) ~ 1 r^L hX vj^ £ 1 , (4) 

requiring the observed radiation to be produced at radii r ra d ^ lO^r^L^i?^ 3 cm 
(cosmological) or r ra d *5 10 7 ^L^r/i 3 cm (halo). 

An additional model requirement is that, for photons observed above pair 
threshold, the "compactness parameter" (or photon-photon optical depth) must be 
r 77 < 1 at the radius r rac [- In the lab-frame (LF) we have r 77 ~ n'^axr/T, where the 
CF photon density in the frame moving with V is = L' / Airr 2 ce 7 = L/47rr 2 ce 7 r. 
Thus 

r 77 ~ La T /(47rrce 7 r 2 ) £ 1 (5) 

for photons observed above the threshold for that T. Otherwise, the above equation 
implies a cutoff of the photon spectrum above e 7 . For the much smaller radii 
r a ~ 10 6 r 6 cm expected as the initial size of the "primary" GRB event from a 
neutron star progenitor, the compactness parameter is large at MeV energies, which 
implies that the initial stages of an impulsive fireball (or for a continuous input, the 
lower portion of the wind) is an optically thick e ± 7 flow, which eventually becomes 
thin at larger radii. 

2. THE GENERIC GRB MODEL 

The observational requirements discussed above provide very strong evidence for 
a relativistically expanding fireball scenario. This has some straigthforward conse- 
quences. The bulk Lorentz factor initially grows linearly with the radius, T oc r, 
until the plasma becomes subrelativistic in its own rest frame. After this T saturates 
to the average value of the initial radiation to rest-mass ratio, r\ = E/Mc 2 = L/Mc 2 , 



r — > 77. The saturation radius r s ~ r/r Q is usually much smaller than the radiation 
radii r ra d required by the (optically thin spectrum) observations. 

One problem faced by simple 1992) fireball models, caused by the satu- 
ration phenomenon, is that beyond the radius r s most of the initial energy has 
been converted to kinetic energy of the baryons, the radiation energy content de- 
creasing adiabaticaly. This would raise enormously the initial energy required to 
explain the observed photon luminosity (Paczyhski , 1990), and in addition it leads 
to photon energy degradation. Another major problem with simple fireball models 
(e.g. Paczyhski , 1986, Shemi and Piran, 1990) is that the only radiation observed 
would be from the fireball photons that escape at the optical thinness radius, and 
these would have a quasi-thermal spectrum. In addition, if the initial event has a 
timescale comparable to a neutron star collapse dynamic time r Q /c ~ 10~ 3 s, the 
observed time over which the fireball becomes thin is also comparable (Goodman, 
1986). Both the spectrum and the timescale would be unacceptable. 

There are, however, at least two ways in which a relativistic fireball can produce 
a longer (10 _1 — 10 3 s) duration, nonthermal 7-ray burst with reasonable energy 
and spectrum. One of them results from the fact that the baryons entrained in 
the ejecta will eventually have to run into an external medium, and there they will 
reconvert their kinetic energy into radiation (e.g. Meszaros , Laguna and Rees, 
1993; Rees and Meszaros , 1992; Katz, 1994). The external medium may either be 
a pre-ejected wind from the progenitor, or the interstellar medium. If its density is 
n Q cm -3 , deceleration occurs at a radius 

r dec ~ lO 17 (£ 51 /n o ) 1/3 (0r/ 2 )- 2 / 3 cm (6) 
and the time-delayed LF duration of the burst is 

tdec ~~ r dec /cT 2 ~ 5 x lO 2 (£ 51 /n o ) 1/3 0- 2/ V 8/3 s. (7) 

Here the total (initial) energy E ~ 10 51 i?5i in a solid angle 9 2 is assumed to be 
released during an intrinsic time shorter than tdec (otherwise, the intrinsic timescale 
t w would be the observed duration). The total photon energy produced in the 
deceleration, for very modest subequipartition magnetic fields which ensure high 
radiative efficiency, is the entire baryon kinetic energy, comparable to the initial 
burst radiative energy. The large radius r rad ~ r^ec ensures optical thiness, and 
the strong deceleration and reverse shocks ensure ideal conditions for relativistic 
particle acceleration leading to synchrotron and inverse Compton (IC) nonthermal 
radiation. 

In addition, shocks may also arise internally in the ejecta itself, before any 
deceleration by the external medium occurs. Such internal shocks could arise for 
various reasons. For instance, the energy or matter input may be time-variable, 
so that higher Y shells overtake lower V shells (e.g Paczyhski and Xu, 1994, Rees 
and Meszaros , 1994). If the energy or matter input at the base of the wind occurs 
during an intrinsic time t w but is modulated on some timescale t v (shorter than t w ) 
with Ar ~ Ar/ ~ T) around an average final Lorentz factor T ~ r], an overtaking 
collision (internal dissipation shock) occurs at 

rdis ~ ct v 7] 2 ~3x lO 14 ^?? 2 . cm - ( 8 ) 
This occurs beyond the wind Thomson photosphere 



r p h ~ M ot I '^TtrripC 2 ''rj 2 ~ 10 12 L 51 r) 2 3 cm 



(9) 



(larger than the saturation radius) for 0.3(L 5 i/t l) ) 1 / 5 ^ r/ 2 ^ 10 2 (L 5 i/t 1 ,) 1 / 4 . Also, 
in general r^is < ^dec- For the magnetic fields turbulently generated in the shocks, or 
for a young pulsar wind, the radiative efficiency of the internal shocks is sufficient to 
radiate an appreciable fraction of the wind bulk kinetic energy. Other mechanisms 
for randomizing the wind energy might be the dissipation of Alfven turbulence 
beyond the photosphere (Thompson, 1994; see also Duncan and Thompson, 1992) 
or the conversion of Poynting flux into photons after the charge density becomes 
insufficient to sustain the frozen-in magnetic field of a young pulsar type wind (Usov, 
1994, 1992). Convective Rayleigh- Taylor instabilities may also arise and become 
nonlinear beyond the saturation radius (Waxmann and Piran, 1994). Below the 
photosphere this would become pressure that reaccelerates the flow, while above 
the photosphere it would be expected to result in freely coasting shells of different 
r\ as discussed above. 

3. SOME OBSERVATIONAL CONSEQUENCES 

What are some of the predictions of the dissipative relativistic fireball scenario? 

One consequence of models based on this scenario is that emission is expected 
at energies other than in gamma-rays. A simultaneous burst of the same duration 
and low but significant fluence is predicted at X-ray and optical energies (Meszaros 
and Rees, 1993b). Breaks in the photon power-law spectrum are also predicted 
(Meszaros , Rees and Papathanassiou, 1994), in qualitative agreement with 7-ray 
spectra (Band, et al. , 1994), and detailed comparison of such calculations with 
observations could provide useful insights into the source physics. For a typical 
MeV fluence F 1 ~ 10~ 6 ergs cm -2 the optical and X-ray fluences predicted for 
cosmological models are compatible with the few detected X-ray flashes as well as 
the lack of widespread X-ray identifications (X-ray paucity), and with the lack (so 
far) of optical detections. Typically , however, they are above the expected HETE 
threshold (Ricker, et al. , 1994) of ~ 10 -10 ergs cm -2 . As an example, for some 
cosmological models in the optical U-band, my ~ 11 — 2.5 log(F M /10 -9 ergs cm _2 ) + 
2.51og(tb/ s), e.g. Meszaros , Rees and Papathanassiou, 1994. Other satellite or 
ground-based observations triggered by BATSE via systems such as BACODYNE 
could, for appropriate sensitivities and pointing times less than the burst duration, 
also detect the GRB at other wavelenghts. Previous attempts at simultaneous 
optical detection (e.g. Vanderspek, et al. , 1994, Krimm, et al. , 1994) did not yet 
have the sensitivity needed for a meaningful comparison, while most other searches 
were not simultaneous. 

A much delayed (days to weeks) radio outburst at the mJy level could also 
become observable when the ejecta has expanded sufficiently for the radio opacity 
to become negligible (Paczyhski and Rhoades, 1993). 

Simultaneous (and delayed) GeV emission is also a fairly straightforward con- 
sequence of this scenario, for electron power-law indices not too steep and moderate 
shock magnetic field strengths (which can be inferred from MeV spectra). A sus- 
tained or delayed GeV emission (as observed, e.g. by Hurley, et al. , 1994) can be 
understood in terms of a wind with internal dissipation and an external decelera- 
tion shock (Meszaros and Rees, 1994). The wind, with an r\ ^ 10 2 and duration 
t w produces via internal shocks an MeV burst extending up to ^ few GeV, and 
also produces a longer event (t ~ td ec , which may be up to hours, see eq[7]) as the 
wind is decelerated in the external medium. At these r\ the external burst is below 
BATSE threshold at MeV, but is prominent at GeV and due to the low compactness 
it extends easily above 30 GeV. An alternative possibility for delayed GeV bursts 
is discussed by Katz, 1994. 



4. DISCUSSION AND PROSPECTS 



A dissipative relativistic fireball scenario is motivated (and practically required) by 
the key observations discussed in §1. It is to a large degree generic, and is expected 
whether the "ultimate source" is, e.g., a young high-field pulsar, a failed supernova 
lb, a neutron star binary merger, a halo neutron star quake, or comets crashing into 
magnetospheres. The basic assumption common to all such sources is that they 
deliver the required energy in a small initial volume (r G 10 6 — 10 7 cm) in a short 
time. Whatever the ultimate source is, it should in any case remain hidden behind 
an opaque e ± 7 veil, and the subsequent evolution of the fireball (during which the 
object manifests itself observationally) is independent of the birth details. This 
may be phrased as a GRB "No-Hair" theorem: the only thing that characterizes 
observationally a GRB is the initial energy, the initial volume and (possibly) the 
initial energy deposition timescale. Knowing the exact nature of the primary GRB 
source would, of course, greatly help in estimating expected rates per galaxy per 
year and details of the spatial distribution. However, this should be (to a good first 
approximation) irrelevant for understanding the physics of the observable intrinsic 
GRB properties. 

A fireball model is, in principle, expected also if GRB arise in the galactic halo. 
In this case the nature, the rate of events and the spatial distribution might be un- 
derstood (see, e.g. Podsiadlowski, Rees and Ruderman, 1995; Colgate and Leonard, 
1994; Wasserman and Salpeter, 1994) if one makes some new assumptions about 
the source physics. In a cosmological setting, on the other hand, there are plausible 
astrophysical sources with relatively uncomplicated physics that could supply the 
observed rate and distribution (e.g. Eichler, et al. , 1989, Narayan, Paczyhski and 
Piran, 1992, Meszaros and Rees, 1992, Woosley, 1993). The dissipative fireball 
scenario described above explains then in a straightforward manner (Meszaros and 
Rees, 1993a, Meszaros , Laguna and Rees, 1993) the GRB energies, the typical 
overall durations and the non-thermal spectra. What is not specifiable in detail in 
such a model is how the energy gets deposited in a low density region to provide 
a high r, which is puzzling either for halo or cosmological sources. However, the 
conclusion that this does occur seems unavoidable, since the observations demand 
such conditions (see §1). 

Several other questions can be addressed within the context of such models. 
One of them is the difference in spectral and temporal properties of halo and cos- 
mological models. The former have smaller total energy and lower magnetic field 
strengths in the wind or at the shocks. If the minimum particle energies accelerated 
in the shocks were very high this could still yield halo GRB spectra which satisfy 
the X-ray paucity and have acceptable MeV breaks, but the time-delayed durations 
are all less than a second (Meszaros and Rees, 1993b). This may be helped if 
the durations are explained with a (so far unspecifiable) internal time t w , but they 
would still tend to overproduce X-rays. In cosmological models, on the other hand, 
either the calculated dynamic time td ec or an internal time t w can give acceptable 
spectra. Details of the spectrum at several wavelenghts, if observed simultaneously, 
could provide discriminants between these cases. 

Another question that may be addressed in these models is the reported bi- 
modal duration distribution (Kouveliotou, et al. , 1993 and these proceedings). One 
possibility is that the durations are given by the deceleration time (7) and bursts 
occur in two main types of external environment. One may speculate, e.g., that the 
progenitor's random spatial velocity causes a fraction of the bursts to occur in the 
galactic halo, while others occur in the disk. The density contrast would be at least 

10~ 3 and the duration difference at least a factor Ahq 1 ^ 3 ~ 10. Alternatively, there 



may be a fraction of bursts for which deceleration occurs in a pre-ejected denser 
wind, while in others the latter is unimportant compared to the ISM. 
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